Construction of Parseval wavelets from redundant filter 

systems 

L. W. Baggett, P. E. T. Jorgensen, K. D. Merrill, and J. A. Packer 

Abstract. We consider wavelets in L 2 (R d ) which have generalized multircs- 
olutions. This means that the initial resolution subspace Vo in L 2 (M. d ) is not 
singly generated. As a result, the representation of the integer lattice Z 
restricted to Vo has a nontrivial multiplicity function. We show how the corre- 
sponding analysis and synthesis for these wavelets can be understood in terms 
of unitary-matrix- valued functions on a torus acting on a certain vector bun- 
dle. Specifically, we show how the wavelet functions on R d can be constructed 
directly from the generalized wavelet filters. 



1. Introduction 

The theory of wavelets is concerned with the Hilbert space I? (R d ) . The prob- 
lem is to find "good" orthonormal bases (ONB), where "good" makes reference to 
several conflicting requirements: 

(1) These bases must be constructed from a small number of model functions, 
called wavelets, and two discrete operations, translation and scaling. In 
this paper, we are concerned with translation by the standard integer 
lattice Z d , and scaling by some prescribed integral matrix A which is 
assumed expansive. 

(2) In passing from function to expansion coefficients, referring to a wavelet 
basis, and back again (this is called analysis/synthesis), the steps must be 
algorithmic, ideally avoiding direct reference to integration over M. d . 

(3) The wavelet functions should have compact support, and should have 
some prescribed number of derivatives. 

The algorithms that have been popular since the mid 1980's are based on 
what is called multiresolution analysis (MRA). This was pioneered by Daubcchies 
[D]. Mallat |Maj and Meyer |Mel [Mey93 , and the idea, while simple, has been 
extremely powerful. The idea itself is much like that of the Gram-Schmidt algorithm 
from Hilbert space, in that it is based on a scale of closed subspaces, resolutions V n , 
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and relative orthogonal complements, detail subspaces W n . The scale of subspaces 
V n play the role of martingales from probability theory. 

Daubechies's book [D] stresses how the requirements Q-© can be met with 
the MRA approach, and all starting with a fixed cleverly chosen function ip in the 
subspace Vb from the resolution. The function tp, the father function, is the solution 
in L 2 to a scaling, or refinement, equation; a solution which results from a cascade 
approximation. The wavelet functions, mother functions, can then be constructed 
from the subspace Wq which is the relative ortho-complement of Vb in V\. One 
drawback of this approach is that if N = |det A\, then N — 1 wavelet generators 
from Wq are needed. In the dyadic case, d = 1, N = 2, that makes one function, 
but in general, N can be large. Now the spaces Vb and Wq are invariant under 
translation by Z d , and there is a corresponding pair of multiplicity functions which 
dictate a minimal choice of generators for Vb and Wo. 

In fact for the general case of L 2 (M d ) and a fixed scaling matrix A, it is possibly 
to get Wq singly generated, i.e., to find a single generator ip. In some cases ip may 
be taken to be the inverse Fourier transform of a subset E of Mr. Such subsets E 
are called wavelet sets, see |DL98t IDLS971 IDLS98) and |B.TMP1LIBMMLIBM| . 
But there are other choices of sets of generating functions ip, with the number 
of generators between 1 and N — 1. What emerges is that these constructions 
force frequency localization, and the compact support in the x-variable is typically 
lost, i.e., we must relax requirement (J3J). As it turns out, the kind of frequency 
localization we obtain is well suited for effective sampling algorithms. A second 
issue enters: The wavelet algorithm may lead to bases which only satisfy a certain 
Parseval property (also called "normalized tight frame"). While we still have the 
resolution structure V n , the number of generators in Vb may increase, but they are 
not directly part of the wavelet basis. This setup is referred to as a generalized 
multiresolution analysis (GMRA); see jB.TMPIL IBMML IBM) . 

The GMRA theory was in fact introduced (in an operator-theoretic/operator- 
algebraic framework) in a pioneering paper by Baggett, Carey, Moran, and Ohring 
BCMO . By now there is a rich journal literature which reflects wavelet construc- 
tions with some degree of multiplicity; see, e.g., |ACML IBLl IDHRSL HTl |RLPS , 
HTCl EH IP^WXI IPTl El IWH [W2 . In addition our multiplic i ty analysis h as 
recently been used in papers by Dutkay and Jorgensen |D.T1L ID.I2L ID.I3L ID.T4j in 
the study of nonlinear dynamical models. 

A recent and interesting paper of Papadakis |P2| offers a global and natural 
approach to generalized multiresolution analyses (GMRAs) based on a geometric 
frame construction which has the advantage of including all GMRA constructions 
in L 2 (R). Some of the differences between the approach in |P2j and the present 
one lie in our use here of operator algebraic tools deriving from the Cuntz algebras 
|Cunj in operator theory. Further, our use of vector bundles in Section 4 offers a 
rather explicit and concrete representation of the matrix functions which in turn 
describe the filters for GMRA normalized tight frame wavelets, alias, for Parseval 
wavelet frames. 

Our paper and other recent papers on wavelet frames may be said to generalize a 
celebrated theorem first proved by Lawton in this journal in the paper |E]. Lawton's 
pioneering result states that a trigonometric low-pass filter which satisfies a certain 
conjugate-mirror filter condition must give rise to a wavelet frame. However, the 
resulting wavelet frame may not be associated with a classical MRA. We generalize 



CONSTRUCTION OF PARSEVAL WAVELETS FROM REDUNDANT FILTER SYSTEMS 3 



Lawton's theorem in several directions to a much broader class of subband wavelet 
filters. Our matrix-valued low-pass filters are still associated with GMRA wavelets. 
These in turn are the most general types of multiresolution structures. 

Finally, we emphasize that the use of wavelet filters derives historically from 
signal processing in communications engineering, see e.g., |Jor03| . 

In the standard case of MRAs, it is well known how the subband filters from 
signal processing allow us to construct the wavelet functions by an elegant algo- 
rithm. The function ip is a solution to a certain refinement equation. The subband 
filters may be thought of as functions on a torus T, frequency response functions. 
But in the case of multiplicity and multiple generators, the corresponding functions 
on T are matrix- valued, and the refinement equation is a matrix equation. 

In this paper we show that the generalized setup admits solutions in L 2 (W. d ) 
starting with this matrix/vector version of the refinement equation. Starting with 
a matrix system of subband filter functions on a torus, we show that our corre- 
sponding wavelet solutions are in L 2 (M d ) , and that they will be Parseval frames 
for the Hilbert space L 2 (R d ) . 

While this is a weakening of the stricter ONB in requirement J2J , the Parseval 
frame property still allows the same recursive analysis/synthesis algorithms popular 
in the MRA case. The now classical method of Mallat and Meyer for constructing 
an orthonormal wavelet in L 2 (M) (relative to translation by integers and dilation by 
2 MaJ IMep proceeds as follows: Let h be a periodic function on K that satisfies 
the "conjugate mirror filter condition" (the so-called Smith-Barnwell condition) 

(1-1) \h(x)\ 2 + \h(x + ^\ 2 =2. 

The function h is the conjugate mirror filter, referred to above, and in our context 
is called a low-pass filter. Consider the infinite product 

p{x) - n -^(j), 

3=1 

and suppose that there exists a nonzero L 2 function <p whose Fourier transform 
cj) coincides with P. Under not too strenuous assumptions on h, this does in fact 
hold. For instance, if we take h to be smooth, and satisfying Cohen's orthogonality 
conditions and the low-pass condition \h(0)\ = ^/2, the set of integer translates of 
<p turn out to be orthonormal functions in L 2 , and in fact is a scaling function 
for a multiresolution analysis {Vj}. 

Given such a low-pass filter function h, there exists an associated periodic 
function g, which also satisfies the filter eauation lI.il and such that h and g satisfy 
the following orthogonality condition: 



(1.2) / l ( x ) 5 ( x ) + / l ( x+ _) 5 ( x+ _) =0 . 

Any such function g, called a high-pass filter, can be obtained from the low-pass fil- 
ter function via the standard technique of constructing a unitary matrix whose first 
row is given by ((h(x)/V2), (h(x + l/2)/V2). Finally, the function ip, defined by 



m = ^(f )?(f ) 
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is an orthonormal wavelet. That is, the collection {ipj,n} = {V% 1 ipi^ 'x — n)}, for j 
and n in Z, forms an orthonormal basis for L 2 (R). 

A famous example of A. Cohen ( jCh) ) shows that eliminating the non- vanishing 
condition can cause the Mallat-Meyer method to go wrong in an interesting but 
not disastrous way. Cohen exhibited a low-pass filter function h, for which the 
infinite product P exists, is the Fourier transform of a nonzero L 2 function <f>, but 
for which the integer translates of 4> are not orthonormal. Further, the translates 
and dilates {ipj >n } of the corresponding function ip, defined just as in the classical 
method, is not an orthonormal basis. So, for this choice of filter h, the Mallat-Meyer 
procedure fails to produce a scaling function, and the resulting function ip is not an 
orthonormal wavelet. Nevertheless, its translates and dilates do form what's called 
a Parseval frame. By definition, this means that for each /6L 2 (R) we have 

ii/ii 2 = EEk/i«i 2 - 

j n 

In P3 and |BJ| , the Mallat-Meyer phenomenon was generalized to incorporate 
the Cohen example in the following way. Suppose h is a low-pass filter for dilation 
by a positive integer N, i.e., satisfies the filter equation 

N-l 



1=0 



and the low-pass condition \h(0)\ — VN, and suppose g\,...,gN-i are corre- 
sponding high-pass filters, i.e., are periodic functions for which the iV functions 
h,gi,..., gN-i satisfy the following orthogonality conditions: 

JV-l 



1=0 



for all 1 < i < N - 1, and 




N-l l 

E St(x + j^)gj{x 

for all i,j between 1 and N — 1. Again, these N — 1 high-pass filters can be con- 
structed from the low-pass filter h by the matrix completion technique. 
Let P be the infinite product 

Then, if the functions h,g\, . . . , gjv-i are Lipschitz continuous, it is shown in BJ 
that there exists a nonzero L 2 function whose Fourier transform coincides with 
P, and the N — 1 functions {tpk}, defined by 

MX) = J=g k ^(±), 



form a Parseval frame multiwavelet. That is, the collection 

{Vi,n,fc = V^MWx - n)}, 
j,n £ Z and 1 < k < N — 1, forms a Parseval frame for L 2 (M.). 
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Just as in the Cohen phenomenon in the Mallat-Meyer constructions, the inte- 
ger translates of the function <f> may or may not be orthonormal, even though their 
closed linear span Vq does form the core subspace of a (generalized) multiresolu- 
tion analysis. Also, the wavelets ipk may or may not have orthonormal translates 
and dilates, and may or may not belong to the subspace V\ G Vq of this associated 
GMRA. 

The arguments in the proof of this result introduce some ideas from operator 
theory, specifically in the form of the Ruelle operator and partial isometries that 
satisfy relations similar to the Cuntz relations |Cun| . In the case of dilation by 
2, W. Lawton [L] used a cascade algorithm in the time domain to independently 
derive the same result (in the special case of trigonometric polynomials, i.e., the 
case corresponding to a single compactly supported scaling function). 

While the constructions of BJJand [L] start with filters associated with a 
classical multiresolution analysis, they build a frame wavelet that may be obtained 
from only a generalized multiresolution analysis. The purpose of the present paper 
is to extend and clarify this kind of result to filters defined in higher dimensional 
and non MRA contexts. That is, we suppose that A is an expansive, integral, 
d x d matrix, and we investigate frame wavelets constructed from filter systems 
associated to generalized multiresolution analyses in L 2 (M. d ) relative to dilation by 
A and translation by lattice points in Z d . 

The theory of generalized conjugate mirror filters relative to a generalized mul- 
tiresolution analysis was first developed in C and BCM . In Section 2 we briefly 
review that subject, and in particular we recall the analogs to Eciuations ll.ll and 
11.21 i.e., generalized filter and orthogonality equations. These generalized equa- 
tions are considerably different from the original ones, because the right-hand side 
now involves an integer-valued multiplicity function on the d-dimensional torus 
T d = [— i, i) d . In general, the low-pass filter becomes a square matrix H = [hij] of 
periodic functions, and the high-pass filter becomes a not necessarily square matrix 
G = [gk.j] of periodic functions. Because the filters in this generalized context are 
matrices of periodic functions, we refer to them as redundant systems of filters. 

By generalizing the arguments in B J] , we are able here to build wavelets from 
generalized filters more simply and with fewer restrictions than in |BCM| . In 
Section 3, under smoothness and low-pass conditions analogous to those in BJ , 
we prove that the infinite matrix product 

oo 

converges pointwise, and the first column $ of P is the Fourier transform of a 
set <j>i, (f>2, ■ ■ ■ of L 2 functions that together generate a generalized multiresolution 
analysis. We then prove that the functions {ipk} given by 

Mx) = ^(^i gftj(A rl W)^ rl (»)) 

form a Parseval wavelet frame for L 2 (R d ). 

In engineering, the various classes of multiresolution analyses are motivated by 
filters from signal and image processing. While these are practical concerns, our 
present aim is theoretical. In most of the engineering MRA-constructions, the mul- 
tiplicity function is equal to 1 on the entire fundamental domain [—1/2, 1/2), so this 
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only produces classical MRA quadrature-mirror filters. One must choose the mul- 
tiplicity function to be mixed, i.e., equal to the characteristic function of delicately 
selected configurations of subsets of [—1/2, 1/2), to get genuinely generalized filters, 
i.e., GMRA filters. If the multiplicity function is not constant, the filters can never 
all be continuous in the frequency domain, which can be a problem for engineers. 
The reason that all the filters for a GMRA cannot be continuous is that they must 
satisfy the filter equations 12.41 and 12.51 both of which involve non-continuous char- 
acteristic functions for non-constant multiplicity functions. This particular feature 
of our theory becomes evident in our example in Section 3, Example 13. 61 Yet, this 
example is surprising in two ways: Firstly, it is related to a non-classical MRA- 
wavelet, the Journe wavelet. Secondly, we construct two low-pass filters which have 
discontinuities in the frequency domain, but the resulting scaling functions are C°° 
in the frequency domain! This prompts the following question: can the scaling func- 
tions and wavelets associated to GMRA's be constructed to have nice properties 
in both the time and the frequency domains? From the viewpoint of engineering , 
this is an important question, and it is addressed in a sequel paper BJMP2 in 
much more detail. We have established this in the affirmative in [B.TMP2] . where 
we construct a frame wavelet that is C' r in the time domain and C°° in the fre- 
quency domain. More importantly, it would be interesting to know whether or not 
scaling functions and frame wavelets that have nice properties in both the time and 
frequency domains can be constructed for every choice of a multiplicity function. 

In BJ] the set of functions h,g\,... , gjv— l is called an M-system. It is shown 
in |B.T| that there is a group that acts freely and transitively on the set of M- 
systems, thus suggesting a natural structure on these systems and therefore on the 
corresponding frame wavelet systems. In Section 4 of this paper, we describe an 
analogous action on the generalized filter systems we have introduced. This time, 
it is a group bundle that acts freely and transitively. 

2. Generalized filters 

We collect here the relevant definitions concerning wavelets and multiresolution 
analyses in L 2 (R d ), relative to translation by lattice points and dilation by an 
expansive integer matrix A; i.e., a matrix each of whose eigenvalues has modulus 
greater than 1. 

Recall that a frame in L 2 (R d ) is a sequence {/„} for which there exist positive 
numbers a and b such that 

o||/ll2<£K/l/»>| 2 <&ll/lla 

n 

for every / £ L 2 (M. d ). It is called a Parseval frame or a normalized tight frame if 
a = b = 1 in the inequalities above. That is, {f n } is a Parseval frame if \\f\\\ = 
EJ</l/n)| 2 forevery/ei 2 (R d ). 

For each z £ TJ 1 , we write t z for the unitary translation operator on L 2 (M. d ) 
given by [t z (/)](£) = f(t + z). Fix an expansive, d x d, integer matrix A and let 
B = A* and N — | det(j4)|. We write S for the unitary dilation operator given by 
[Hf )](.*-) = V^V/(A(f)). For each element lj G T d , there exist exactly N distinct 
points ( G T d such that = cu, where a denotes the endomorphism on T d in- 
duced by the action of B on M. d . We denote these points ( in a Borel manner as 
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u>o,uji,. . . , uin-1- Note also that, because A is an expansive matrix, the endomor- 
phism a is ergodic. 

Definition 2.1. A multiresolution analysis (MRA) in L 2 (R d ), relative to the 
group {r z } of translation operators and the dilation operator 6, is a sequence 
{V^}^oo of closed subspaces of L 2 (M. d ) for which: 

(1) VjCV j+1 . 

(2) V s+ i=S{V s ). 

(3) [JVj is dense in L 2 (R d ), and tlVj = {0}. 

(4) There exists an element <f> G Vq whose translates {t z (</>)} form an or- 
thonormal basis for Vq . Such an element (f> is called a scaling function for 
the multiresolution analysis. 

A generalized multiresolution analysis (GMRA) is a sequence {Vj} of closed 
subspaces of L 2 (R d ) that satisfies conditions (1), (2), and (3) above, but satisfies 
the weaker fourth condition 

4'. Vq is invariant under all the translation operators r z . 

In both cases, the subspace Vq is called the core subspace of the GMRA. 

The theory of GMRA's is developed in B CMO] . jBMMj. and [BM] (see 
also |Plj h In particular, it is shown there that by Stone's Theorem for unitary 
representations of the group Z d , or more generally by spectral multiplicity theory for 
a set of commuting unitary operators, the invariance under translation of the core 
subspace Vq implies the existence of a multiplicity function m : T d >— > {0, 1, 2, • • • oo}. 
The multiplicity function m counts the number of times each exponential function is 
represented as a subrepresentation of translation by Z d on Vq. We will assume in this 
paper that m is bounded, with maximum value c. We let S t = {lu eT d : m{u) > i}, 
and recall from BM that there exist generalized scaling functions {4>i}i<i< c such 
that the collection {T z (4>i)} for z £ Jj d and 1 < i < c forms a Parseval frame for 
Vb, and such that J2zez d |<Mw + 2:)| 2 = xs^w). Note that these results translate to 
the classical conditions in an MRA, which is the special case of a GMRA in which 
the multiplicity function is identically 1. 

If Vq is the core subspace of a generalized multiresolution analysis {Vj}, then 
the subspace Wo = Vj. fl Vq- of a GMRA {Vj} also is invariant under all the trans- 
lation operators t z . Hence, again by Stone's theorem, there exists a complementary 
multiplicity function rh : T d {0, 1, 2, • • • oo} that characterizes the representation 
of Z d on Wq. As a direct result of the fact that Vi = Vq ffi Wq, the multiplicity 
functions m and fh must satisfy the following consistency equation (see BMM ): 

N-l 

(2.1) m(u>) + m{u) = m{ui). 

1=0 

By the consistency equation, the assumption that m is bounded implies that m is 
bounded as well. We write c for the maximum value of fh, and Sk for {ui £ T d : 
rh(uj) > k}. 

Generalized multiresolution analyses are useful because of their relationship to 
wavelets. In particular, in jBMMj it is shown that every orthonormal multiwavelet 
is associated to a GMRA with m =the number of wavelets. We recall the relevant 
definitions here: 



s 
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Definition 2.2. An orthonormal wavelet, or more generally an orthonormal 
multiwavelet, in L 2 (R d ), relative to the group {t z } of translation operators and 
the dilation operator S, is a finite set ipi, . . . ,ip m of functions in L 2 (R d ) for which 
the collection {4>j,z,k] = '(Tz(V'fc))} forms an orthonormal basis for L 2 (R d ). A 
set of functions if>\, . . . ,if> m forms a frame multiwavelet if the collection {ipj,z.k} 
forms a frame for L 2 (M. d ), and the set is called a Parseval frame multiwavelet if the 
collection {ipj tZ ,k} forms a Parseval frame for L 2 (R d ). 

We now are ready to develop the definition of generalized filters. The classical 
filter equation fEauation ll.lJ) . which is the basis for the Mallat-Meyer construction 
of scaling functions and wavelets, emerges naturally from the study of classical 
multiresolution analyses. For, if is a function, for which the closure of the linear 
span of the translates of <p is invariant under the dilation f(x) — > /(x/2), e.g., if <f> 
is a scaling function for an MRA, then it follows from elementary Fourier analysis 
that <p satisfies the following classical refinement equation: 

(2.2) \/2<j){2x) = h(x)$(x), 

where h is a periodic function that satisfies the filter equation 1.1. As described in 
the introduction, the Mallat-Meyer construction |Maj . [Me] reverses this procedure, 
by beginning with a sufficiently well behaved function h that satisfies filter equation 
1.1, and producing a corresponding scaling function and multiresolution analysis 
by iterating eauation l2.2l to get an infinite product expression for 4>. 

In an analogous way, the theory of generalized filters emerges naturally from 
the study of generalized multiresolution analyses (see (£] and |BC M ). Because 
V- 1 = <5 -1 (Vo) is contained in Vq, it follows that, for each i, there exists a sequence 
{hij} of functions on the torus T d such that 

(2.3) y/Nfc(B(x)) = J2 h *A x )^ x )- 

3 

As is shown in |BCM| , these generalized filters satisfy the generalized filter equa- 
tion: 

N-l 

(2.4) Yl h iAui) h i',M) = k^ N xsM, 

j 2=0 

and also have hij supported on Sj. In analogy with the classical case, we make the 
following definition. 

Definition 2.3. A generalized low-pass filter relative to a GMRA with multi- 
plicity functions m and m is a matrix of functions H = [hij] on T d (or equivalently 
7L d periodic functions on R d ), with hij supported on Sj (or the periodization of Sj), 
that satisfy both the generalized filter equation (|2.4|l and the generalized low-pass 
condition hi,j(G) = Si^Sj^y/N 

Just as in the classical case, we can sometimes reverse the procedure of pro- 
ducing filters from wavelets. In particular, generalized low-pass filters give rise 
to generalized high-pass filters. The relevant result is contained in the following 
theorem, again from BCM . 

Theorem 2.4. Let H = [hi.j] be a generalized low-pass filter relative to a 
GMRA with multiplicity functions m and fh. Then there exists a matrix G = [gk,j] 
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of functions on T d satisfying 

iV-l 

(2-5) 9k,j(ui)gk',j(ui) = 5k : k'Nx-g k {uj), 

j 1=0 

and 

N-l 

(2-6) J2 J2 h ^ ; (^)9^M) = 

j 1=0 

for all i and k. The matrix of functions G — [gk,j] * s called a generalized high-pass 
filter 

Under conditions which allow the production of generalized scaling functions 
from generalized low-pass filters, these high-pass filters can be used to build frame 
multiwavelets. Some narrow conditions that allow this are described in (BCM| . In 
Section 3, we show that more general conditions allow the same construction. 

There is one final property of a generalized low-pass and high-pass filter that 
we will need in Section 3, which was not presented in BCM . 

Theorem 2.5. Let H and G be a generalized low-pass and high-pass filter 
relative to multiplicity and complementary multiplicity functions m and rh. Assume 
that the maximum value of m is c and the maximum value of rh is c. Then 

c c 

(2.7) ^2hi d (uJi)hi,j'(uJi>) + ^2g k ,j(uJi)gk,j'(uJi') = Sj^Si^NxSji^i)- 

i=l k=l 

PROOF. For each oj 6 T d , we define a matrix K (u>) having c + c rows and 
C x N columns as follows: We index the c x N columns of K(u>) by pairs (j, I), 
where 1 < j < c and < I < N — 1. Then the entry k^MuS) is defined to be 
(l/^/N)hij(uji) if 1 < i < c, and ki^j^(uo) = (l/\/iV)gi_ C j(a>f) if c < i < c. We see 
directly from Equation H2.4(l that, for 1 < i < c, the ith row of the matrix K{u>) 
contains a nonzero entry if and only if u> £ Si, i.e., if and only if i < m(w). And, 
from Equation 12. 5f) . for c < i < c + c, the ith row of K(lu) contains a nonzero 
entry if and only if lj £ Sti— C \, i.e., if and only if i — c < m(cu). Therefore, there are 
exactly m{uj) + rn(w) nonzero rows in K{lj). 

Next, we note that the column indexed by the pair (j, /) has a nonzero entry 
only when some hi^(uji) or gij(ui) is nonzero. That is, the (j, /)th column has a 
nonzero entry only when loi G Sj, i.e., only when j < m{uji). So, the maximum 
number of columns having a nonzero entry in them is equal to X^o 1 TO ( tL 'i)j which, 
by the Consistency Equation l]2.ip . equals m(uj) + rh{oj), and this is exactly the 
number of rows of K(u) that have a nonzero entry. 

Therefore, the set of nonzero entries in the matrix K(u) are contained in a 
square submatrix L(uj) of size (m(w) + fh(uj)) x (m(w) + m(w)). 

Finally, from Equations (|2.5(l . l|2.4l and l|2.6|l . we see that the rows of this square 
matrix L(uj) are orthonormal. Hence, the columns of L(u>) are also orthonormal, 
and this implies the orthogonality equations of the theorem. □ 

3. Construction of frame wavelets from generalized filters 

We are now ready to use the generalized filters from Section 2 to extend the 
construction procedure for frame wavelets described in |B.T| . Just as in the classical 
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case, the first step of the construction is to take an infinite product of dilations of 
the low-pass filters. 

Proposition 3.1. Let H = [/ijj] be a generalized low-pass filter as in defini- 
tion \2.!A Assume that the components of H are Lipschitz continuous functions in a 
neighborhood o/O. 



(1) The infinite product 



q=l V 



converges almost everywhere on R d , and the entries P^j of P are square- 
integrable functions on M. d , with Pjj = for j > 1. 
(2) For 1 < i < c, let & be the inverse Fourier transform of the function P»,i. 
Then the <pi 's are continuous at 0, satisfy 0^1(0) — 5i i, and also satisfy 
the following generalized refinement equation. 



1 c 

<f>i(B(x)) = -=^2h^{x)cj)j(x). 



3=1 

Proof. Throughout the proof, we will use the following result from linear 
algebra. If C is a c x c matrix all of whose eigenvalues have modulus less than 1, 
then J^fcLo ll^' fc ( a; )ll converges, for every x E R c . This result follows by using the 
Jordan canonical form for the matrix C to show that ^ ||C fc (ir)|| is dominated by 
a convergent geometric series. 

We prove the convergence of the infinite product first. If (3 is the Lipschitz 
constant, note that by the low-pass condition we have < P\\x\\ for {i,j} ^ 

{1,1}, and ||hi,i(a:) - VN\\ < /3\\x\\, both for ||x|| sufficiently small. Write P k for 
the partial product n g =i "^^(^ _9 ( 2 '))' wnere as before TV = |det(A)|. We first 
show by induction that for each fixed x, there is a bound K such that \P^j(x)\ < K 
for all 1 < i,j < c, and k > 1. To see this, write 



hi 



1=1 v 



E^r x (*) 



1=2 



-4n s " fc wii 



N 



for k sufficiently large (where in the first term we use the fact that the orthogonality 
conditions give |/ii,i| < yN). Using induction on k, and the linear algebra result 
mentioned above, we get the bound we seek. 

Now, using this bound, we see that for fixed x and for j > 2, \P k ,^(x)\ < 

YTi=i Pif 1 ^) ~jjy\\B~~ k ( x )\\ ^ as fc ^ oo. For the remaining case of j = 1, we 
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have 



\PtM-Pk\x)\ 

1=1 v^V 

= Pti 1 ^) (j=h 1 , 1 {B- k x)-l) +j^P^\x) (-L hlil (B- k 3 



< 



Ei^*)i-4ii* 



N 



< -Z=cK\\B- k x\\. 



We then see that {P k \(x)} is Cauchy and thus convergent, again by using the linear 
algebra result mentioned at the beginning of the proof. 

To complete the proof of (1), it remains to show that the limiting functions 
Pi i are in L 2 (R d ). To do this, we will first use induction to prove that 



j=l J B k (Q) 



P k Ax)\ 2 dx < 1. 



the Pjj are periodic modulo 



We note that because the hkj are periodic modulo Z' 
B k (Z d ), and thus the domain of integration can be replaced by any set that is 
congruent to B k {Q) modulo B k (Z d ). To select the replacement set, we first choose 
coset representatives ui\ = 0, u 2l ■ ■ ■ ,wjv of Z d /B(Z d ). Since 



N 



J (Q + L0 n ) = B(Q) modulo B(Z d ), 



n=l 
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we can take as our domain of integration the set B^ 1 (\_\ n=1 (Q + u n )) Using this, 
we have 

E / \P&(x)\ 2 dx 

~[ J B k (Q) 

= e/ N {iiPti\x)^hAB- k X )\ x 



\m — 1 V 

N k ~ 1 ( £ P*r 1 {B k - 1 x)P£- 1 (B k - 1 x) 



l x) dx 



—=hij(B 1 x) ■ h m i(B~ 1 x)dx 

s/N J WN jy ' 



r c N i 

J( 3 2,m j=l n=l ViV 



h m< j(B 1 x — B 1 ui n )dx. 



We note now that modulo Z d , the set {-B l u>i, B 1 ui2, ■ ■ ■ ,B 1 wn} parame- 
terize the distinct N elements of Q = [— |, i) d that map to (0, 0, • • • , 0) under the 
cndomorphism a : T d ^ T d , where here we are identifying T d and \) d . Using 
the orthogonality relations satisfied by the hij, the last equation simplifies to 



x)| dx 

(Q) ' 



N"- 1 ( Y,Pti\B k ^x)P^\B k ^x) XSl {x)dx < E/ \ P i,T\ 

= E / i^r'wi 2 ^- 

It follows that 

E/ i^wi^^i 

~J J B k (Q) 

for all G N, and since U^L -B fc (Q) = R d , by Fatou's Lemma we have 

E / \PiAx)\*dx<l. 

In particular we get Pij E L 2 (R d ), which completes the proof of part (1). 

The refinement equation in part (2) is immediate. It also follows from the proof 
above that the infinite product P converges uniformly on neighborhoods of 0, and 
thus that the ifii's are continuous at 0. Finally <pi(0) — <5,.i is a consequence of the 
low-pass condition. □ 

We will now use the results of Proposition 13.11 to build a frame wavelet. Wc 
begin by generalizing the computational ideas in |BJ| . Define two Hilbert spaces H 
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and H by H = 0° =1 L 2 (S l ) and H = 0£ =1 L 2 (S k ); and two operators S H ■U^H 
and S G : H -> H by: 

^(/JlH^'H/faH) 

and 

[s G (/)K^) = gV)/(«M), 

where / = /j G 7i, if* and G* denote the transposes of the matrix functions H 
and G respectively, and as above a denotes the map on the torus T d induced by 
the action of the transpose B of A acting on R d . Note that fj(a(ui)) and fk(a(oj)) 
are necessarily only defined when a(ui) belongs to Sj for the first case and Sk for 
the second. We define these functions to be outside of these domains. We record 
here the formulas for the adjoints of the two operators: 



(3.1) 
and 



1 



[SM/)]M = ^E ff ^)/(^), 



z=o 



N-l 



z=o 



S* H 



It will also be convenient to have explicit formulas for the powers of both Sh and 

n-l 

[SS(/)]H = HHV( W ))/KH, 



fc=0 



and 
(3.2) 



N" 



1=0 fe=n-l 

where the w;.„'s are the N n points £ € T d for which w = a n (£). 

We also include here an estimate involving these operators that we will need 
later, 



(3.3) 



c c N n — 1 

EE E 

i=l j'=l Z=0 



n-l 



I]H*(a k K n )) 



Lfc=0 



< cTV" 



which we prove using induction. The case n = 1, 

c c jV-1 

]T£]TfeMI 2 <civ, 

i=l j = l 1=0 

follows from equation II2.4JI together with the fact that J2j=i XSj (w) = m(uj) < c. 
For the induction step, note that the w;,n+i can be labeled in such a way that 
a n (uJsN+q,n+i) = Wq, so that uj S N+q,n+i = {u q )s,N- Thus, writing I = sN + q and 



14 L. W. BAGGETT, P. E. T. JORGENSEN, K. D. MERRILL, AND J. A. PACKER 



using the Cauchy-Schwarz inequality, we obtain: 

2 



c c AT™ + 1 -1 

EE E 

i=l 3=1 



1=0 
c N 



k=0 



< 



< 



EE E 

i=l j=l i=0 
c c AT"-1JV-1 

EE E E 

i=l j=l s=0 g=0 



IJ ff*(a fc ( WJin+1 )) 




n-1 



k=0 




E 



Lemma 3.2. The operators Sh and Sg satisfy the following relations: 

(1) S h Sh = I, the identity operator on TL. 

(2) SqSg = I, the identity operator on TL. 

(3) S* h Sg = and S g Sh — 0, the operators on the appropriate Hilbert 
spaces. 

(4) ShS^ + SgS g = I, the identity operator on TL. 

Proof. These are direct consequences of the orthogonality relations satisfied 
by generalized filter systems relative to m and rh. For example, to prove (1), for 
/ G H, we write 



S* H S H f(u) 



1 ^ 



1=0 



N-l 



i££/uv)//'uv-|/u-) 

1=0 



where the last equality follows from 12.41 The other parts of the lemma follow 
similarly from l2~5l 12.61 and 12 . 71 respectively. □ 

We note that both Sh and Sq are partial isometries, and that the relations in 
the lemma are similar to the famous relations defining the Cuntz algebra Oi- We 
now use these operators to decompose the Hilbert space TL in a convenient way. 

Lemma 3.3. Let Rq denote the range of Sg, and let R n — S r ^(Ro). Then TL = 
©^Lo^"- Moreover, if z is any element in Z d , and fk, z is the element ofTL whose 
kth component is the exponential function e 27Tl ( u \ z ) and whose other components 
are 0, then the collection {S H (Sa(fk,z))}> for k running from 1 to c and each z 
running through Z d , forms a Parseval frame for the subspace R n . Consequently, the 
collection {S H (SG(fk,z))}, 1 < fc < c, z g Z d , and < n < oo, forms a Parseval 
frame for TL. 

Proof. That the subspaces {R n } are orthogonal follows directly from the 
relations in Lemma 3.2. That the elements {S H (S g (fk,z))} form a Parseval frame 
for R n follows from the fact that Sh and Sg are partial isometries, together with 
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the fact that the elements {fk,z}, as k runs from 1 to c and z runs over Z d , form a 
Parseval frame for 7i. 

Write Rao for the orthogonal complement in H of the direct sum (J) R n . We 
must show that i?oo = {0}. Note that Roc is invariant under both Sh and S H , and 
therefore that the restriction of both these operators to Roo are unitary operators. 

By way of contradiction, suppose that /o is a unit vector in Roc. For each 
natural number n, write /„ = S H n (fo). Note that f n +m = S* H m {f n )- We make two 
observations about /„. First of all, for almost all u> € T d , we have 

||/„HH 2 = ||sT(/oM)ll 2 

c c N n —i r o 
= ^EEE II H{a\^ n )) [/„]>,,„) 

p=l 1=0 lk=n— 1 



< 



iV 



i=l 

^ c c N"-l 

J" EE E 



2,p 



c JV™-1 



n ^(« fe (^,«)) 



n— 1 



c AT — 1 



E i[/o] P '(^',«) 



p' = l l'=0 



< C X 



j^E E l[/o] P '(^',n)| 2 , 



where the last inequality follows from the transpose of equation EH By the point- 
wise ergodic theorem this then implies that 



limsup||/ n H|| 2 <c f \\f (w)\\ 
Jo 



doj = c. 



Secondly, since Sh is unitary on Roc, we have 



N"-l 



— H/oKr. 



N n -1 



1=0 



± \\s H s H n (fo)M\\ 



1=0 
N n -1 c 



E 



< 



]\[n EZ EZ 

i=o i=i 

^ c JV"-1 c 

jyn EZ EZ EZ 

i=l i=0 p=l 
2 



n-1 



fc=0 
n-l 



[S£ B (/o)]pM 

2 

xEIWp'HI 5 



< c||/ n (w) 

implying (again by the pointwisc ergodic theorem) that 



liminf ||/„(o;)|| 2 > 



1 



Consequently, by Egorov's Theorem, for any e > 0, there exists an Mo and a 
set E C T d such that the measure of the complement of E is less than e, and 
i — e < ||/ n (tj)|| 2 < c + e for all w S -E and all n > Mo,. Before we apply this 
theorem, we will establish some other estimates needed in our choice of e. 
First, we pick an integer Kq > 3 log 2 c + 9, so that we then have 



(3.4) 



K +l 3 

N >32c*. 



(This follows since N = | det A\ > 2.) 
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Next, we choose a 5 > so that 
(3.5) 5 < 

and 
(3.6) 



1 



4c K °+ 3 VN 



N(l-S) 



K + l 



Sc K ° 



> 



N 
2 



Note the second condition is possible since the function on the left hand side ap- 
proaches y/N as 6 goes to 0. 

From the low-pass condition and the requirement of Lipschitz near 0, we know 
that for any 5 > there exists a neig hborhood U of e T d such that \h itj (uj)\ < S 
for all u) G U and all pairs ^ (1, 1), and \hi^(ui)\ > VN(1 - S) for all uj G U. 
Let U be the neighborhood corresponding to our choice of 5 satisfying (1) and (2) 
above. By continuity of a, there further exists a neighborhood V C U such that 
for every wGVwe have a k (uo) G U for all < k < 2{Kq + 1). Hence, if P(u>) is the 
matrix given by P(ui) = l\l=o Ht {a k {uj)), then for f (1, 1), < I < K + 1 
and lj &V, 

\?iAo?{u))\<6{c^ K °. 

(There are c K ° summands, each having Kq factors, and in each, one factor is 
bounded by S, and the other factors fbv l2.4fl are bounded by Vn.) For (i, j) = (1, 1), 
we have 



(VN(1 - S)) Ko+1 - 5{cVN) Ko < \P 1A (a l (iu))\ < VN K "^ L + S(cVn) 

(Again, there are c K ° summands, each having Kq + 1 factors. One of these sum- 
mands is bounded below by (y/N(l — 5)) Ko+1 and above by y/N 0+1 , and the other 

(c K ° — 1 summands are bounded by Ss/N °.) 

Now choose an e smaller than the measure of V, and also smaller than l/(4c), 
so that we are assured that the corresponding set E will satisfy X(E n V) > 0. In 
fact, we may even assume that the set of u> E V for which a k (oj) € E n U for all 
< k < 2(Ko + 1) has positive measure. 

Fix an u G E n V for which u k (uj) eEnUiorO<k< 2(K + 1). If M 
is the natural number corresponding to this choice of e as above, and n > Mq, 
we have the following estimates on the components of f n {w). First, for i ^ 1 and 
< I < K Q + 1, 



K 



l/n>'M)l 



[5f 0+1 (/„ 



+k +i 



^P lJ (a'(c))[/„ +i , 0+ i],(^+^+ 1 ( W )) 
i=i 



< 



E 

j=i 



■Pij(o'M) x |[/„ + K +i]j(a 



i+ATo + l 



M)| 



< c<5(cVA0 °%/c+~e 



so that by condition l|3.5|l on 6, 
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El[/n]i(a'M)| 2 < 4c£V(*°+ 2 )iV*° < I. 
Therefore, because u> € E, for all n > Nq we must have, for < / < l^o + 1, 

|[/n]i(a'(u,))| 2 >^-e>l 

On the other hand, 

c 

[/ n ]r( W ) = ^P lj (^)[/„ + K 0+ i],(a^ +1 (c)) 

C 

- P 1 , 1 (o;)[/„ + K 0+1 ] 1 (a^ +1 M) +^P 1J M[/ n+ K 0+ i] J -(a^ +1 (a;)), 

3=2 

implying that 



|[/„ + k 0+ i]i(^+V))I < 



|[/n]l(«)| +Ei=2 |^jH| >< |[/n+K +l]i(a ifo+1 H)| 



< 



< 



< 



< 



IA.iMI 

+ c<5(cv / AQ J< " v / c+7 

(VJV(1 - - 6(cy/N) K ° 

Vc+e 



■if o+l 



N 
8v^ 



■if o + l 



1 

4? 



where in the third step we use conditions (|3.5|l and (|3.6|) to simplify the numerator 
and denominator respectively, and in the final step, we use condition (|3.4|) . Hence, 
the point a Ko+1 (uj) satisfies 

I < \if n+Ko+ i}i(a Ko+i n)\ 2 < ±- c , 



Ac 



which is a contradiction. 



□ 



We now state our main result. 



Theorem 3.4. Let fa be defined from the infinite product of the low-pass filter 
system H as in Provosition \3.1\ and let G = [gk,j] be the corresponding high-pass 
filter system. Then, if we define a function ipk € L 2 (U. d ), for 1 < k < c , by 

V 3=1 

the collection {ipk} is a Parseval frame wavelet for L 2 (M. d ) relative to dilation by A 
and translation by lattice points z G Z d . Further, if Vq is the closed linear span of 
the translates of the fa 's by lattice elements z € "L d , then {Vj} = {<^(Vo)}jgz is a 
generalized multiresolution analysis for L 2 (M. d ). 
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PROOF. We prove first that the {V'fc} form a frame wavelet. For convenience 
in what follows, we introduce the following notation. For each 1 < k < c and each 
z € Z d , write r y k ' z for the element Sa{fk,z) of Rq C H, where as in the previous 
lemma, fk, z is the element of Ji whose fcth component is the exponential function 
e 2TTi(ui\z) an j w h ose ther components are 0. Note that ^ k ' z = 0j =1 lj' Z , where 

2iTi(a(cj)\z) 



7,- '*(<*>) = 9k,j(u)e 
It then follows from Lemma 3.3 that for each F € Tt we have 

oo c 



ii f ii 2 = EEE i(^i^(7 M ))i 2 . 

n=0 fe=l zgZ d 

For each integer n, each 1 < k < c , and each z e Z d , define the function tp n k,. 



by 



^„, M (i) = \/]v n V fe (^ n (i) + ^, 

and note that the Fourier transform of "4> n ,k,z is given by 

c 

^T z ( x ) = AT ^ e 2 ^ B ~ n (*> I z > ]T ff fc j (B" n - 1 (a:) )^(B" n - 1 (*) ) . 

j=i 

We wish to prove that the collection {V'n.fe,^} is a Parseval frame for L 2 (R d ). 

Now, let / be an element of L 2 (M. d ) whose Fourier transform / has compact 
support. For a fixed integer J > 0, we have 

E EEk/i^m>i 2 

n=-oo fe=l z gZ d 

= E EEk/i^)I 2 

n=— oo fe=l z gZ d 



- E EE 

n=— oo fe=l zeZ"" 



E EE^ 



n— — oo fc— 1 . 



/ /(x)JV-r- 1 V]e 2wi < B " n ( a! )l z >flfc J -(B- n - 1 (a;))0j(B- n - 1 (a;))da; 
/ /(B J+1 (x))e 2 ^< s -" +J+1 Wk> x 



-M+2J+1 



= EEE^ l+1+J 

n=0 fc=l z gZ d 



2 5fc j(B-"+-'(x))0 J -(B-"+-'(x))tfe 



/ /> J+, (x)) e MB»+ 1 Wk)V ij (B»( I ))^(B»( I )) ( li 

j=1 



= EEE^ 1+J 

n=0 fe=l zgZ d 



/ /(B J+1 (i)e«( B " +1 WWV 5tj (B»(i)) 



Yl H(BP{x)<S>(x) 

_p—n—l 



dx 
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EE E ^ 1+J | /.EE J{B J+ \x + O)^ + 0T[Ht{BP{x)—z(B n (x))dx{< 

n=0 k=l Z £l d Td j=l C&L d P=0 



n=0 k=l z gZ d 



IF 



J||2 



H ' 



where F J — (J) F J j is the element of TC given by 



F J ^) = VN 1+J E f(B 1+J (x + C))4>j(x + 0- 



Hence. 



E EEk/i*. 

n= — oo k=l z£E d 



k.z) 



7=1 ^ 



E 



j=i j B1+ , (T£!) cezd 

E / XB 1 + -'(T d )( ; c) 

3 . =1 yR- 



J2 K x + B 1+J (C))0i(5- 1 - J (x) + C)| 2 ^ 



Now, because / has compact support, and the matrix B is expansive, there exists 
a J' such that the support of / is contained in B 1+J (Q). There must also exist, by 
the compactness of B 1+J (Q), a Jo such that B 1+J (Q) contains B 1+J (Q) for all 
J > Jq. Now, for J > Jo, the product Xb 1 + j (Q)( x )f( x + B 1+J (Q) is nonzero only 
if x = B 1+J {y) for some y G Q, and also x + B 1+J {() = B 1+J (y + () £ B 1+J {Q). 
Consequently, for any J > Jo, and any ( ^ 0, we must have Xb 1 + j (Q) ( x )fi x + 
B 1+J (0) = for a11 x + Hence, for J > J we have 

E EE k/i^,m)i 2 = E I \?{ x )?\TAB- l - J { x ))\ 2 dx, 

n=-oo fc=l z gZ<i 1 = 1 "' Rd 



so that, by the Dominated Convergence Theorem, we obtain 

oo 



E EEk/i^ 

n=-oo fc=l Z £Z d 



|/(x)| 2 y^|^(0)| 2 dx 



l/ll 



This demonstrates the Parseval frame property for elements / whose Fourier trans- 
forms have compact support. For general / S L 2 (R d ), the Parseval equality follows 
from the density of these functions. 

It remains to show that {Vj} form a GMRA. Properties (2) and (4') of Definition 
12.11 are direct consequences of the definition of the Vj, and property (1) follows 
immediately from Proposition 13.11 (2). The density required in Property (3) is a 
consequence of the fact that the {ip n ,k,z} form a Parseval frame for L 2 (M. d ). 
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To finish the proof, we now show that DVj = {0}. Write {$} for a sequence 
of elements whose translates form a Parseval frame for Vo (such a sequence must 
exist since Vo is closed under translation), and Pj for the orthogonal projection 
operator onto the subspace V}. To prove that flVj = {0}, it will suffice to show 
that lirn^oo ||P_j(/)|| = for each / € L 2 (M. d ). By a standard approximation 
argument, it will suffice to show this holds on a dense subset of L 2 (R d ). Thus, let 
/ be a Schwartz function for which / vanishes in some neighborhood Nf of 0, and 
write C f for the (finite) number £ fe |/* /* (fc)l- Such f's are dense in L 2 (R d ). The 
Poisson Summation Formula holds for such an /, and we will use it in the following 
form: 

E \K B ~ j tt + 0)| 2 = E / * f*( Aj (k))e- 2mm) ■ 

l k 

Now, for each £ S T d , let be the largest number for which B~ J + I) E Nf 
for all || 1 1| < Because A (and thus B = A r ) is expansive, we must have that 

tends to infinity for almost every £. Finally, we use the fact that the function 

i i I 

is integrable on the cube. Hence, we have 



\P-Af)f 



^EE 

i z 

EE 



-2™<B'(£)|z) 



-T 

No 



-E 

No ^ 



Yd 



E m+i)f(B- j (t+i)) 

\m\>h(z) 

^ jvjE/ E m+Q\ 2 E \f(B- j (t;+i))\ 2 dz 



/ 

Jt 

J 

Jt 



< 



E E m+w 

i \i\>ij(0 

E E mz+W 



/EE $fc+oi 
<c/ EE i^ + oi 2 de. 



l 



Ei/*-r(^'( fc ))i 



which approaches as j goes to infinity by the Dominated Convergence Theorem, 
the integrand here being bounded by to. □ 
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Theorem 13 .41 shows that the functions ipk are in the set Vi, so that the Parseval 
frame wavelet {ipk} that we have constructed is obtained from the GMRA {Vj}, 
in the sense defined by Zalik Zl . ft is an open question whether the wavelet is 
associated with the GMRA {Vj} in the sense that Vj is the closure of the span of 
{S l (T z (tpk))}i<j (see e.g. |Bw| ) . Moreover, the multiplicity function of the GMRA 
produced by the theorem may not coincide with the multiplicity function used in 
the construction of the filters, as the first of the examples below (Example 13.5(1 
shows. 

Example 3.5. Let d = 1, A — 2, and the multiplicity functions m and fh both 
be identically 1, so that we start in the classical setting. For our filters, we take 



X-i,i)+X±[J,|) 



h = V2 

and 

9 = V2 \x ±hi) + X±[hi) 

Note that h and g satisfy the definitions of generalized low- and high-pass filters, 
and are Lipschitz in a neighborhood of 0. The infinite product P — YfjLi "^5 M^r) 
equals X[— 1 \)- The integer translates of the function <j> whose Fourier transform is 
this infinite product P, are not orthonormal, and do not determine the core sub- 
space of any classical multiresolution analysis. On the other hand, as predicted 
by Theorcm l3.4l the standard construction in this case does produce a generalized 
multiresolution analysis {Vj}, with multiplicity function m = i). As guaran- 
teed by Theorem 13.41 the construction also produces a Parseval frame wavelet tp, 
here given by 

#z)=X[-i,-i)+X[ii)- 

In this case the function (j> is easily constructed out of negative dilates of if), so the 
wavelet is necessarily associated with as well as obtained by the GMRA. 

We end this section with an another example for dilation by 2 in L 2 (IR). 
This one begins with filters for the multiplicity function of the well-known Journe 
wavelet, whose Fourier transform is the characteristic function of the set 

r 16 ^ r 1 2, r 2 1, . 16, 

[-y,-2)U[--,--)U[-,-]U[2, y ). 

We use Theorem l3.4l to build a Parseval wavelet with a C°° Fourier transform that 
is associated with the Journe multiplicity function. 

Example 3.6. The Journe multiplicity function m takes on the values 0,1, 
and 2, with Si = [-5, -f)U[-|, |)U[|, ±) and S 2 = j)- The complementary 
multiplicity function fh{x) = 1, since the Journe wavelet is a single orthonormal 
wavelet. 

To build the filters, we let po be a C°° classical (MRA) low-pass filter for 
dilation by 2 (i.e., that satisfies the classical filter equation II. 1|) . that in addition 
satisfies po(%) = for x e ±(\ - e, + e) U (f - e, | + e). Note that byOl we then 
have po(x) — \/2 for x £ ±(= — e, ^ + e) U (— — e, ^ + e). Let pi be the standard 

choice of C°° high-pass filter associated to po, given by pi{x) = e 27Ttx po(x + |). 
Then it is easily checked that the following functions satisfy 12.41 12.51 and 12.61 and 
thus arc generalized low- and high-pass filters by our definitions in Section 2. 
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hi,i(x) 




9i(x) 
92(x) = 



\ Po (x) ze[-f,f) 

1 otherwise 

Po(x+i) x&[-\,\) 
otherwise 

= (V2 xe±[f,i) 

1 otherwise 
h 2 , 2 (x) = 

Jpi(z) xe[-f,|) 
| otherwise 



\pi{x + \) xt[=±,\) 
I otherwise 
Now we check that the resulting wavelet does in fact have a C°° transform. 
We first fix an x and show that <j>i and <j> 2 are C°° in a neighborhood of x. 
Recall that these functions form the first column of the infinite product matrix 

Yl'jLi ^7f ^/j 1 '*! * j ^ 1 ' 2 q 23 ^ . Because fi2,i is on (— |, |), all but a finite num- 
ber of the matrix factors are upper triangular. Thus, each of the entries in the 
first column of the infinite product contains only a finite number of terms. Since 
hi.i = \pl in a neighborhood of 0, each of the terms has only a finite number of 
factors not equal to 1. Thus, it will suffice to show that each of the factors in each 
of the terms is C°° . 

By construction, we have that the hij are all C°° everywhere except for h\ t \ 
at n ± |, ft.1,2 at n ± i, and h 2 A at n ± | for n e Z. We will show that whenever 
one of these discontinuities occurs as a factor in the infinite product, it is cancelled 
by a following factor that is in a neighborhood of the point of discontinuity. Note 
that the product of two adjacent factors in the infinite matrix product has the form 

) + /ii, 2 (2/)/i2,i(!) /ii,i(y)/ii,2(|) N 
/i2,i(y)/ii,i(f) h 2 . 1 { y )h 1 . 2 {^) / 

Thus any term in the infinite product that contains a factor of h 2 ,\(n ± |) must 
also contain a factor of one of the forms /ii,i(n ± |), h\^{n ± ^), hi, 2 (n ± 
h\. 2 {n ± |). The first three possibilities are in a neighborhood of the point in 
question, so if we have a discontinuous factor of /i2,i, it is either cancelled out by 
a factor, or we also have a factor of h\.\(n ± |) with a smaller n. Similarly, any 
term in the infinite product that contains a factor of /ii i i(n± |) must also contain a 
factor of one of the forms /ii j i(n± yj), /ii j i(n± ? ), hi t2 (n± or hi, 2 {n± j). The 
first three of these possibilities are in a neighborhood of the points in question, so 
any discontinuous factor of h\ t \ is either cancelled out by a factor, or is followed 
by a factor of h\^ 2 (n ± i) with an equal or smaller n. Finally, any term in the 
infinite product that contains a factor of h\^ 2 (n ± ^) must also contain a factor of 
either the form h 2 .i(n± j^) or the form h 2 .i(n± |). The first of these possibilities 
is in again in a neighborhood of the point in question; the second possibility 
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throws us back into the first type of discontinuity we considered above, but with 
a smaller n. We can repeat the above sequence of arguments if necessary, noting 
that each succeeding factor is evaluated at a point half the distance from the origin 
as its predecessor, so that the chain above must eventually end with a factor of 0. 

This argument shows that cf>i and 4>2 are C°°. To see that ip is C°° as well, it 
suffices to note that g\ and g 2 are. 

Remark 3.7. With some more work it is possible to slightly improve the ar- 
guments above and to show that in fact the Fourier transform vanishes rapidly at 
infinity. 

4. A Generalized Loop Groupoid Action on the Bundle of Generalized 

Filter Systems 

Let {hij}i< i t j < c and {<?fc,j}i< k <5, i< j <c be generalized low-pass and high- 
pass filter functions defined as in Section 2. Since ®f =1 L 2 (Si) = L 2 (\_\^ =1 Si), 
we can suppress the second index of the filter functions and view generalized filter 
functions as a vector (c + c-tuple) of functions: 

c 

(fti.fta, ••• h c: g 1: g 2 ,--- ,9b) G [L 2 (\JS t )} c + d . 

i=l 

Further, we note that for any fixed u> £ T d , all the information in the output of the 
vector of functions 

(hi,h 2 , ■ ■ ■ h c ,gx,g 2 , • • • ,35) is actually in c m ^ u ^ +0i( - a ^\ where as in Section 
2, a is the endomorphism of T d defined by the matrix B — A f , since by the orthog- 
onality relations, hi(u>) = if i > m(a(u))) and gk(x) = if k > m(a(u))). 

We want to generalize the discussion given by Bratteli and Jorgensen in [B.Tlj , 
and construct a loop groupoid which acts on the class of filter systems with bounded 
multiplicity functions corresponding to a fixed dilation matrix A acting on R" . We 
will also impose the condition the filter systems satisfy some "initial conditions" 
that are in some sense canonical. 

We first remind readers of the notion of a vector bundle over a topological space 
X; more details can be found in |FD| . 

Definition 4.1. Let X be a topological space. A (finite dimensional) vector 
bundle over the space X, denoted by (E, p, X), is a topological space E, together 
with a continuous open surjection p : E — > X, and operations and norms making 
each fiber E x = p~ l {X) into a (finite dimensional) vector space, which in addition 
satisfies the following conditions: 

(i) |/H || is continuous from E to R, 

(ii) The operation + is continuous as a function from {(?/, z) G E x E : p(y) = 
p{z)} to E. 

(iii) For each A G C, the map y 1— > A ■ y is continuous from E to E. 

(iv) If x € X and {yi} is any net of elements of E such that \\yi\\ — > and 
p(y%) — > X in X, then y t — > G E x in E. 

A Borel map s : X — > E is called a Borel cross-section if po s(x) — x, Vx G X. 

We review the notion of essentially bounded multiplicity functions m associated 
to a dilation matrix A that can give rise to GMRA's, as described in |BCM| and 
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|BM] , We first recall that to must satisfy the consistency inequality 

JV-l 

(4.1) m(w) < Y m (^i), 

1=0 

where we recall from Section 2 that {loi : < I < N — 1} is the set of N 
distinct preimages of ui E T d under the endomorphism a given by a(ui) — B(ui), 
for uj E T d = R. d /Z d and B = A 1 . This inequality leads to the consistency equality 
Eauation l2.1l discussed in Section 2. Moreover, we can enumerate the tui as follows. 
Enumerate a set of coset representatives {£o = 0, £i, • • • £/v-i} for B~ 1 (Z d )/Z d 
including E B-\Z d ). For each I, < I < N - 1} let Q be the image of under 
the natural projection from R d onto T = M. d /Z d . Note that the Q are distinct 
elements of T d and that a(&) = 0, < I < N - 1. Find a Borel cross-section 
a : T" = R d /Z d -» T™ = M d /B^ 1 {Z d ) =i [R d /Z d ]/[B- 1 {Z d )/Z% with a(0) = 0. 
This Borel map satisfies ao a(ui) = u> for all ui E T d . Then define 

wi = ct(w)G, < I < JV - 1; 

one easily verifies that the {ui : < I < N — 1} are the N distinct preimages of 
uj E T d under a. 

Let A = U^L B k (S! + Z d ). Recall from [EM] and BRS that in order for m 
to be the multiplicity function for a GMRA, A must satisfy 

(4.2) Y XA(w + n) > m(w). 
Finally, A should satisfy 

(4.3) U peZ £P(A) =R n . 

If conditions 14.11 14.21 14.31 are satisfied we say that m is an (essentially) bounded 
multiplicity function with respect to the dilation matrix A. Given such an to, we 
construct the conjugate multiplicity function m by defining 

N-l 

(4.4) m{u>) = Y m M - m M, w € T d ; 

1=0 

by definition, m and fh will satisfy Eauation l2.ll 

Given an (essentially) bounded multiplicity function m on T d , let c = ess sup to, 
and c = ess sup to. We recall from Baggett, Courter and Merrill that to such an 
to we can explicitly construct a canonical system of low-pass filters, or "generalized 
conjugate mirror filters" {h% j : 1 < i,j < c}, by using the method of Theorem 
3.6 of [BCM]. Given this system {htfj : 1 < i, j c}, Theorem 2.5 item (1) of [BCM] 
gives us a way of explicitly constructing an associated system of high-pass filters, 
or "complementary conjugate mirror filters", {g% • : 1 < k < c, 1 < j < c}. We 
call this family of filters {h^ : 1 < i,j c}U {g% ■ : 1 < k < c, 1 < j < c} the 
canonical filter system associated to the multiplicity function m. 

Now let T. } = {lo e T d : m(a(w)) + fh{a(J)) = j}, < j < c + c. Set 
T itj = S l nT 3 , < j < c + c; then each T itj is Borel and Si = \Jj=l T i,3- 
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Definition 4.2. Fix a bounded multiplicity function to associated to the di- 
lation matrix A. Let E m be the Borel space given by 

c c+c 

/•,, U IJP^ xCi ]- 

i=l j=0 

Let (E m , p, Uj =1 Sj) be the Borel vector bundle where the map p : £" m — > Ui=i Si 
is defined by p(x,v) = x, (x,v) £ Ti_j x O. By definition, (E m , p, S*) is 

a vector bundle over |J^ =1 Si whose fiber over a; S Si is a complex vector space 
of dimension m(a(w)) + m(a(u>)). An M-system associated to the multiplic- 
ity function to is a Borel cross-section M : Uj=i Sj — > £"m of this bundle 
whose values, (Mi(u>), M2(cj), • • • M m ( Q ( (iJ )) + m(a(w)))i are the output of a vector 
of generalized low- and high-pass filters. That is, for fixed u> £ U^=i Sj> the first 
m(a(u>)) components of the vector M(u>) correspond exactly to the generalized 
low-pass filter function values hi(oj), h,2(u>), • • • , /i m (( a ( w ))(w), and the last m(a(w)) 
components of M(u>) correspond to the generalized high-pass filter function val- 
ues gi(u>),g2{u)), ■ ■ ■ , <7m(a(a;)) where the filters involved satisfy the low-pass 
and Lipschitz conditions defined in Section 2, and where in addition the filters 
{hi}\j{gk} satisfy the initial conditions /i$(C/) = h^(Q), 1 < « < c, < I < N — 1, 
and <7fc(0) = fffc(O)) 1<^<S, < Z < JV — 1, where the /if and the g% are the 
canonical filters associated to m in the previous paragraph. Denote by Ai m the set 
of M-systems associated to the bounded multiplicity function m. We remark that 
M m can be given the structure of a topological space if its elements are viewed as 
elements of the Hilbert space ®= =0 L 2 (T hj ) C j . 

Remark 4.3. The definition above is essentially the same as Definition 4.2 of 
BJ MPl] except for the additional assumption that the M-system have the ap- 
propriate canonical values at the preimages of under the endomorphism a, which 
was missing from BJMPl . Note all the information about the generalized filters 
{hij} and {gkj} is encoded in the M-system. In particular, for any fixed multi- 
plicity function m, such that both m and rh are constant in a neighborhood of the 
origin, we have a one-to-one correspondence between M-systems and collections 
of generalized filter functions as defined in Section 2 which satisfy the canonical 
initial conditions. These generalized filters will in turn give rise to an orthonormal 
frame wavelet family and its associated GMRA {Vj} by Theorem 13 . 41 from Section 
3. However, as shown by one of the examples given in Section 2, the multiplicity 
function m! for the GMRA {Vj} need not be equal to m. So, there is no correspon- 
dence in general, between the class of M-systems associated to a given multiplicity 
function m and the class of GMRA's associated to to. 

To develop the loop group action on these M-systems, we first define an en- 
domorphism n Q : |Ji=iSi -> T d by IL a (oj) = a(uj). Each u> e T d has 
^^q 1 m(u)i) = m(aj) + to(c<j) preimages in Si, where the {uji : < I < N— 1} 
are the N distinct preimages of ui in T d under the endomorphism a, parametrized 
as discussed earlier. For convenience of notation, we label these preimage maps 
rnj<\, where rnj\(uj) = u)\ G Sj C Ll° =1 Sj for 1 < j < m{u>i). (Note that this 
range on j, as I varies from to N — 1, gives all the preimages, since if j > m(u>i), 
by definition loi is not an element of Sj.) For each fixed w, we give the pairs 
the lexicographical order, so that (h,ji) < (fejia) if h < h or if l\ = Z 2 and j\ < 
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We thus implicitly define a 1-1 map A w taking the pairs onto the integers from 
1 to m(oj) + rh(uj). 

We now construct a unitary group bundle (F, q, T d ) as follows. For each 
j € {l,--- ,c + c, } } let Zj = {uj £ T d : m(uj) + m(uj) = j}. Let 

C+C 

F m = \_\[Z j x U(j,C% 

3=0 

where U(J, C) is viewed as a the topological group of unitary j x j matrices: it 
inherits its topology from being a subset of C*-algebra of complex j x j matrices 
given the operator norm. Define q : F m — > T d by g(w, T) = w, for (w,T) € F m , 
and note that q : F m — > T d is a continuous open surjection, and the fiber g -1 ^) 
of the bundle consists of the group of complex unitary matrices U (m(uj) + fh(uj) , C). 
Borel cross sections to this group bundle consist of Borel maps L : T d -> F m such 
that g o L(u>) — uj. Note also that (F m , q, T d ) is a subset of a Borel vector bundle 
over T d in the sense of Definition 14. II whose fibers over particular values of uj E T d 
are made up of algebras of finite-dimensional matrices of varying dimension. Wc 
denote the set of sections of this bundle by T m (F m , q). Note T m (F m , q) is a group 
under pointwise operations on T d , where the identity element of the group is given 
by that section whose value at uj is equal to /d m ( u )+m(u)' 

Also, it is possible to define a groupoid corresponding to a dilation matrix A as 
follows. Let VLa consist of the set of all bounded multiplicity functions correspond- 
ing to the dilation matrix A, that is, 

flA = {miT^NU {0} : m is essentially bounded and satisties ll.ll4.2l4.3ji -. 

where two multiplicity functions for A are identified if they are equal almost every- 
where on T d . There are a variety of topologies we can put on f2^i; for the moment the 
Hilbert space topology of L 2 functions on T d will do. We call SIa the multiplicity 
function space associated to A. Now set 

Ca = |_J T m (F m ,q). 

Then Ca is a groupoid, with range map r equal to the source map s defined from 
Ca to fl a by r(j) = s(j) = m for 7 6 T m (F m ,q). Hence two elements 71 and 72 
of Ca 

can be multiplied together if and only if ^1, ^2 £ r m (F m , q) for a fixed 
bounded multiplicity function m. 

We are now state a theorem about M-systems that can be derived in a fairly 
straightforward fashion from the orthogonality relations. This theorem is a gener- 
alization of Theorem 4.3 of BJMP1 , and thus we merely sketch the proof and 
refer the reader to BJMP1 for details. 

Theorem 4.4. Let T m (F m , q) be the group of cross sections of the group bun- 
dle associated to a fixed bounded multiplicity function m defined above. Let M : 
Sj — > E be an M-system associated to m. Then uj h- > (ii,A„(i,j) where 

is an element ofT m (F m ,q). 
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Proof. As noted above, 1 < X u (l, j) < m(uS) + rh(uS), so for each lu € T, the 
matrix (L^x^i ) {to)) is a square matrix of the correct dimension. 

We shall show that for all uj 6 T d , the rows of (Li x u iij)(u>)) are orthonormal, 
by means of the orthogonality relations for generalized filter functions given in 
Equation |^7| in the statement of Theorem 12. 51 

Write Li for the ith row of (L^wy)^))- If 1 w(w), 

m(a;)+m(aj) 



< Li, Li> > = ^i,Au,(i,j')( £J )-^«»^«('»j)( w ) 
A w (i,j)=l 



^Mi(r (i (w)) J ^-Mi> (r { i j) (a;)) 



At_i m(wj) 

X! E \ V -M.,,.,M-My y 
1=0 3=1 

N-lm(u)i) 

y - 

i=o j=i 

j c 2V-1 



- E 



hij(ijj{)hi',j(uji) 



N ^ ^ hi :j (cji)hi> ;j(ui) 

3=1 1=0 

(as hi t j(u>i) = for j > n(u>i), since ^ Sj in that case) 

= (by Equation (2.2)) N^^xsM = *M' 

(we note that (w) = 1 since we have i < m(o>) and for those values of i, to G S 1 ; 
by definition of m(u>)). 

The cases m(ui) m(tj) + m(w) and 1 < i < m(uj) < i' < 

m(w) + m(w) follow from similar arguments using Equations 2.4 and 2.5. Thus we 
have that in all cases, the rows of (Li t \ u n fi((j)) are orthonormal, and we have the 
desired unitary matrix. □ 

Remark 4.5. We note that for uj e T d , the (m(uj) + m(uj)) x (m(uj) + fh(uj)) 
matrix L(u>) is exactly the submatrix L(u>) of the matrix K(u>) defined in Theorem 
2.5. Thus in some sense the proof given above is redundant. We have added in 
the extra detail because we want to exercise care in enumerating the rows and 
columns of L(ui) for future use. We also note that the "initial condition" on M- 
systems given in Definition 14.21 can be rephrased as follows. Let L c iX j)( u )) t> e 
the element of T m (F m , q) corresponding to the "canonical" filter system associated 
to the multiplicity function m, and let Lq be be the [m(0) +m(0)] x [m(0) + fh(0)] 
unitary matrix L c iX n j) (0) • Then for any M system M associated to m, we have 

y^M,-(r(jj-)(0)) = Lq. This follows from the initial conditions satisfied by any 
M system outlined in Definition 14.21 

The results of Theorem 14.41 imply that the columns of (-ki,A„(Z,j)( w )) are or- 
thonormal as well, and allow us to obtain the following corollary, which is a gener- 
alization of Corollary 4.4 of BJMP1 . 

COROLLARY 4.6. Letm andrh be bounded multiplicity and conjugate multiplic- 
ity functions associated to a dilation matrix A that are constant in a neighborhood 
of the origin, with related sequences of sets {Si | 1 < i < c} and {Sk | 1 < k < c}. 
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Suppose {hij}i< ij < c and {<?fc,j}i< k <e, i< j <c are generalized low-pass and high- 
pass filter functions associated to the multiplicity function m. Then for all ui G T d , 
and for all I, V G {0, 1, • • • , N — 1} and j, j' G {1, • • • , c)}, we have 

C+C 

£ 77 M i( r ('j)( w )) M i( r C',J')M) 



= ^j^h ld {uji)h ld '(uji) + ^2jjgk,j(ui)9k,j>(ui) = Sjj'Sij', 

i=l fe=l 

where here the Mi correspond to the generalized filter functions hi y j and gkj as 
given in Definition \4-ty 

PROOF. This result follows fairly directly from Theorems 2.5 and 14.41 and we 
omit the details of the proof. □ 

We are ready to define the generalized loop group and its associated action on 
the set of M-systems associated to a multiplicity function to. 

Definition 4.7. The loop group associated to the multiplicity function to is 
defined to be the subgroup Loop m (F m , q) of the group of Borel sections T rn (F m , q) 
whose elements K satisfy L(0) = W m (o)+m(o)> anc ^ ^j,A„(i,0 are LiPi in a neigh- 
borhood of the origin. 

We now prove that the generalized loop group Loop m (F, q) acts freely and 
transitively on the set of M-systems associated to the multiplicity function to. The 
following theorem generalizes Theorem 4.7 of BJMP1 . 

Theorem 4.8. There is a free and transitive action of Loop m (F, q) on the set of 
M-systems associated to an essentially bounded multiplicity function m such that to 
is constant in neighborhoods o/0/, < I < N — 1, where the set {0; : < I < N — 1} 
consists of the preimages of under the endomorphism a : T d — * T d This action 
is given by 

L-M{u) = L(n a H)[(Af 1 H,M 2 (c J ),--- ,M m(QM)+ ^ (QM) (c J ))] t . 

Proof. We first prove the transitivity. Suppose we are given two different 
M-systems, labeled M = (Mi) and M = (Mi). Define an element L of the group 
bundle associated to to, that is, an element of T m (F, q), where L(u>) has dimension 
m(uj) + m(u>), as follows: 

^ m(cj)~\-rh(uj) 

(4.5) L it i>(uj) = — Yl M i'(r(Zj)M)M(r (i j)(w)). 

Since the M-systems M and M have the same values at and at the preimages Q of 
under the automorphism a, one easily verifies that L(0) — W m (o)+A(o) ■ Also, one 
sees by inspection that if Lm and L-^ are the elements of T m (F m , q) corresponding 
to the M-systems M and M as in Theorem 14.41 then the proposed element L of 
Loop m (F, q) given in Equation IOJ is exactly L(u) = L-jj(w)[Ljvf (w)]*. 
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In addition, as shown by Bratteli and Jorgensen in the classical case, we have 

m(oc(u)))-\-fh,(a(uj)) 

[L-MUuj) = J2 Li,i'(naM)Mi,(w) 



i' = l 

m(a(uj))+m(a(ctj)) / m(a(cj))-\-rh(a(uj)) 



E 4 E M,(r (!j) (a( W )))M,(r (lj . ) ( tt (,))) M* 



TV 

A Q( c)(/j) = l 

m(a(a;))+m(a(u;)) / m(a(a;))+m(a(ct;)) 

53 M(r (u) (a(u,))) ^ M v (r m (a{u)))M if {w) 

Aa(u,)(i,j)=l V = 1 

= Mi(«), 

where the last equality follows from the orthogonality of the columns of M as 
established in Corollarv l2.7l so that the inside sum is except for the single values 
of I and j where rij(a(u>)) = lu. 

To prove that the action is free, suppose M = (Mi) is a M-system associated 
to m and L 6 Loop TO (F, q) satisfies 

L(Il a (u)))[(M 1 (uj),M 2 (uj),--- ,M m(QM)+ ^ (aM) (u;))] t = 

[(M 1 (u),M 2 (lj), ■ ■ ■ ,M ro(a(w))+jJi(aM) (w))]*. 

Using Theorem l4.4l for each ui £ T d we define a (m(a(u)))+m(a(uj))) x (m(a(ui)) + 
rh(a(ui))) unitary matrix M. by 

M ^ a ^)(U3)(u) = J^Mi(r (lJ) (a(u))), 

1 < i < m(a(uj)) + m(a(u)), < I < N - 1, 1< j < m(a(uj)i)). 

We then see that L(a(uj))M(uj) = M(u>) for all u e T d . By unitarity of M(oj), 
this shows that that L(a(uj)) is the (m(a(u>)) + m(a(u>))) x (m(a(aj)) + m(ai(u;))) 
identity matrix for all ui G T d , which implies that L is the identity element of 
Loop m (F, q), as desired. □ 

Similarly, we can define the generalized loop groupoid assigned to a dilation 
matrix A. 

Definition 4.9. Let A be a d x d integer dilation matrix. Let Qa be the 
multiplicity function space associated to A. The loop groupoid associated to A 
consists of the set 

Coop A = \_\ Lo °Pm( F m,q), 

where r — s : Coop a — > is defined by r(L) — s(L) — m for L S Loop m (F m , g). 

We can also define a fiber bundle of M-systems over f2,4 as follows: 

Definition 4.10. Let VLa be the multiplicity function space associated to A. 
Define M.a = U me n A {m} x M m , where M m consists of the set of M-systems as- 
sociated to a fixed multiplicity function m, where two M-systems are identified if 
they are equal almost everywhere. Give elements of M m the Hilbert space topol- 
ogy mentioned in Definition 14.21 The topology on is obtained from viewing an 
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element (m,A4) 6 Ma as a tuple (m,M.) of functions defined on the Cartesian 
product space T d x U^JT*]*, taking on values in [N U {0}] © ? 2 (N), where 

(m,M)(uj,z) = (m(^),M 1 (z),M 2 (z), - ■■ M m ( a ( z )) +7h ( a ( z) )(z), 0,0,0, - ■ ■), 

and where M.j{z) = Mj(z) if z G |_li=i ^« an< l 1 < J < m(a(z)) + m(a(z)), and 
= if z Ll^_i for {Si} the standard subsets of T d associated to the 
multiplicity function m. Then Ma is a topological space, if we view elements of Ma 
as elements in the Hilbert space L 2 (T d ) © L 2 (\J^ 1 [T d ] i ) <g> [Z 2 (N)]. Define a map 
II : — > by H((m, A4)) = m. Then II is a continuous surjection, since it is 
a restriction of the Hilbert space projection from L 2 (T d ) © £ 2 (|_|ti \^ d ]i) ® [^(N)] 
onto L 2 (T d ) © [{0}] = L 2 (T d ), and n -1 (m) = M m , so that MU is a fiber bundle 
over n^, called the fiber bundle of A/-systems associated to A. 

Using these definitions, we obtain the following corollary to Theorem 14.81 

Corollary 4.11. Let A be a d x d integer dilation matrix, and let Coop a and 
M.a the loop groupoid and fiber bundle of M -systems associated to A, respectively. 
Then there is a groupoid action of Coop a on Ma, where L € Coop a is allowed to act 
on M £ Ma if and only if r(L) = H(M) = m. This action is fberwise transitive. 

Proof. This is just a restatement of part of Theorem 14.81 □ 

We now modify Example 4.5 of BJMP1 , in order to show how we can use the 
loop group action to transform the canonical filter functions for the Journe wavelet 
into the filter functions discussed in Example 13. 61 The construction of the canonical 
filter functions for the Journe wavelet was first done in the thesis of J. Courter 

Example 4.12. The Journe wavelet in the frequency domain is the character- 
istic function of the set 

hy,-2)U [--,--) U [-,-] U [2, y). 

Here the multiplicity function m takes on the values 0, 1, and 2, and m(x) = 1, since 
the Journe wavelet is a single orthonormal wavelet. If we identify T with 

we can write Si = -|)U [-f , f )U [f , i), S 2 = [-±,±), and Si = [-£, i]. The 
canonical generalized filter functions then are: 

hi,i(x) = \/2X[-f,-J)u[-i ,i)u[i 

gf(x) = n/2x[_i _!)□[!,!) (a;); 

7*1,2 (*) = 0. 
h c 2a {x) = 0, 

g%(x) = V^Xl-ii)^)- 
One calculates that the "initial condition" matrix Lq discussed in Remark 14.51 
corresponding to this canonical filter system is the 3x3 matrix 



L c = ( 1 

1 
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We now construct an element of the loop group L p such that if Mj is the M-system 
corresponding to the above output of generalized low- and high-pass filters, the M- 
system L p ■ Mj has as its output the filter functions corresponding to those given in 
the example in Section 3. Consider the decomposition of the disjoint union Si |J S 2 
(identified with a subset of [— 5, |) |J[— |, |)) given by 



T1.1 


= ± 




j£) (Here m(x) > 1, m(2x) + rh(2x) = 


1-) 




Tl,2 


= ± 


V 

14 


i) U ±[jj, |) (Here m(x) > 1, m(2x) 
U ±[f , |) (Here m(x) > 1, m{2x) 


■f fh{2x) 


= 2 


71,3 




1 

14 


4- rh(2x) 


= 3 


72,2 


= ± 


14 


A) (Here m(x) = 2, m(2x) + m(2x) = 


2.) 




72,3 


= h 


1 

14 


jj) (Here m(x) = 2, m(2x) + rh(2x) - 


= 3.) 





Then as in the general case, Sj — |_|^ =0 7},fc, j = 1, 2. 



We now describe the M-system Mj : Ui=i &i ~ ¥ Em associated to the filter 
functions above. We describe Mj as a cross-section separately on both of the 
(disjoint) sets Si and S 2 . 
On Si , Mj is given by: 



[Mj{x)\ 



= < 



[V2], if seTi.1, 

[^2,0], if xe±[ii)U±[i,f)CT li2l 

[0,V2], if x e ±[£, i) C Ti, 2 , 

[0,v^,0], if a; e Ti, 3 \T 2 ,3, 

[A 0,0], ifx€[Ti )3 nTa 1 3]. 



On S 2 , Mj is given by: 

[MjOb)] 



[0,V2], ifxeT 2 , 2 , 
[0,0, V2], if.TGT 2 , 3 . 



Consider the decomposition of the circle T (identified with [— i, A)) given by 

Pi = [-i, i) (Here m{x) = 2, m(|) = 2, m(^±i) = 1.) 

P 2 = ±[i, |) (Here m(x) = 1, m(f ) = 2, m(*£l) = 0.) 

P 3 = ±[|, f ), (Here m(x) = 0, m(|) = 1, m(^±i) = 0.) 

P 4 = ±[f , |) (Here m(x) = 1, m(f ) = 1, m(£±I) = 1.) 

We note the following, which will be useful in our calculations: 27\ 1 = P 3 , 
2(±[A,i)) = 2(±[I,f)) = P 4 , 2Pi, 3 = Pi, 2(±[^,I)) - 2T 2 , 2 = P 2 , and 
2T 2 , 3 = Pi. 

Now as in Example l3.5l let po be any classical (MRA) low-pass filter for dilation 

by 2, (i.e., one that satisfies the classical filter equation ll.lfl . that also satisfies 

p (x) = for x £ ±(A - e, + e) U (f - e, f + e). Note that by Equation HT1 we 

then have p Q {x) = ^2 for x e ±(f - e, £ + e) U (-i - e, £ + e). Let Po = ^| 

be the normalization of po, so that Po(0) = 1. Let pi be the standard choice of 

associated high-pass filter, so that pi{x) — e 2mx p (x + |), and let p\ = ^= be 

the normalization of p\ . Again, by Equation 11.11 p\ must have modulus 1 for x € 
±(I-e,A +e)u (3_ e) 4 + e) . 
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The associated element of the loop group L p that we want to choose is: 

if x € Pi, 



L p (x) 



P~o(f) pB(2±i) 




1 







pl^ii) 



P~o(f) Po(^) 
Pi(f) Pi(^) 
Pi(f), 



pw: 



if a; G P 2 , 
if x 6 P3, 
if x G P 4 . 



On Si, our new M-system [M pJ (x)} = [L p (2x)[M J (x)] t ] t is the M-system M Pr 



defined by 



[pi (a;)], if x G Ti,i, 
[po(a;),pi(x)], if x G ±[£, i) U ±[^, |) C T 1)2 , 

[ Po (x), P i(x)}, if a; G±[|,f) CTi, 2> 

[Mp.jKx) = <( [0,\/2,0], if x G ±[§, |) = Ti, 3 \T 2 ,3, 

[po 0*0,0, Pi 0*0] > if x G [-i, i) = Ti :3 n T 2 ,3, 

\po(x+ |),pi(a;+ ifxGT 2! 2, 

, [Po(z + i),0,pi(a;+ §)], ifxGT 2 , 3 . 

We thus obtain the generalized filter functions coming from M p .j : 
fri.ifc) = Po(aOX[-f,f)0*0, 



^2,1 0*0 

- 

h p 12 {x) 
S 2 0*0 



Pi(x)X[-§ i)u[i ,1)0*0, 
= Po(a: + -)x [ _i^ ) (i), 



/if, (a;) 



0. 



Pi(a: + 

Note that these are exactly the filter functions obtained in Example 13. 61 

One could no doubt adapt the above example to obtain generalized filters sim- 
ilar to those given in Examples 4.2 and 4.3 of BCM . 

Acknowledgments. The authors gratefully acknowledge helpful conversa- 
tions with Astrid An Huef and Iain Raeburn. 
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